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The six elastic constants (and six elastic compliances) of rutile were deterniiiied in the 
kilocycle per second frequency range by a resonance method. The slanchud (h'viations 
range from 0.2 percent for sn to 4.3 percent for S13. 



1. Introduction 

The electrical properties of rutile make it a prom- 
ising material for solid state electronic devices and 
a large literature on such properties exists and has 
been surveyed [/, 2, 3].^ (\)mparatively little infor- 
mation on th(^ mechanical j)rop(Mlies of rutile was 
available until recently. The linear compressibilities 
had been measured by Bridgman [4], and a set of 
elastic constants was calculated from spectroscopic 
data by Dayal and Ap])alaTnirasimliam [5]. No 
further work on elastic properties was apparently 
done until the recent independent and nearly simul- 
taneous measurements of three groups. Vick, Hol- 
lander, and Brown [6], determined four of the six 
elastic constants of rutile directly from pulse ve- 
locity measurements and calculated the other two 
by combining Bridgman 's linear compressibilities 
with their own data. Verma [7] published a set of 
the elastic constants calculated from his pulse ve- 
locity measurements. The present writers reported 
[8] a set of values determined by least squares fitting 
of the theoretical equations for the orientation de- 
pendence for elastic moduli to a set of experimental 
values determined on 10 single-crystal rutile rods by 
resonance in the audiofrequency range. Appreciable 
discrepancies existed for some of the constants, es- 
pecially between those of Verma and those of the 
other investigators. Vick and Hollander [9] checked 
their measurements and published a refined set of 
values which are little different from their first set. 
The most serious discrepancies were removed when 
it was discovered that the direction which Verma 
took to be [100] was in fact [110]. It is an interest- 
ing property of the 4/mmm crystalographic point 
group, to Avliich rutile belongs, that the x axis may 
be tidven as either [100] or [110] and a self consistent 
matrix of elastic consttmts results which has the 
same form (but three of the six constants will have 
different numerical values). Ttiere is thus no in- 
ternal inconsistency in Verma's work but his results 
should be expressed in terms of [100] as the x axis 
to agree with the usual convention and to permit 
direct com])arison with other values. This calcula- 
tion has been done by Birch [JO] ami the resulting 
values agree fairly well with the results of other in- 

1 Italicized fisuros in brackets indicate the literature references at the end of 
ithis paper. 



vestigators but not as well as would b(^ (wpected from 
the accuracy of methods of measurement which were 
used. The present writers wished to refine their 
data; their initial results were based upon measure- 
ments on 10 rods whose rod axes wqtq all near the 
[001] axis. A comparatively large uncertamty in Cn 
and C06 resulted. This situation was improved when 
SIX additional rods were eventually obtained and the 
present results agree well with the values calculated 
by Birch in the sense that there is now no statistically 
significant difference for any of the c^^ Tlie stand- 
ard deviations are, however, appreciably largcn' tlian 
those obtained by the authors on corundum [11] 
using the same resonance method. The most im- 
portant source of this variability is believed to be 
caused by the presence of a profusion of small angle 
boundaries whose existence is shown by back-reflec- 
tion Lane ])at terns. 

The present paper gives i\w method of applying 
the resonance teclmique to tetragonal crystals 
because it apparently has not been previously 
described for this crystal system. '^Fhc^ general 
method has been described previously and reference 
[11] may be considted for general background. 

2. Description of Specimens 

All specimens were sjnitlietic ^'single crystals'^ 
grown by the Verneuil flanie-fusion technique using 
an oxygen-hydrogen flame. Such specimens in the 
as grown state characteristically contain two types 
of defect: (1) they are oxygen deficient and (2) they 
have many small angle boundaries (as much as 2° 
misorientation across a single boundary). The first 
type of defect is easily removed by heating in oxj^gen 
at temperatures around 1,000 °C. Parentheticall}^ 
we note that the difficult problems of obtaining exact 
stoichiometry and of measuring small deviations from 
stoichiometry do not concern us here. The elastic 
constants do undoubtedly depend to some extent on 
iJie degree of reduction but the effect is very small 
as noted by Vick and Hollander [9]. We have not 
carried out a systematic study of tlie effect of reduc- 
tion on all of the Ctj but observations of the very 
small changes produced in several heavily reduced 
(jet black) rods indicate that tlie effect which the 
possible small remaining oxygen deficiency in our 
nominally stoichiometric rods may have on the Cij 
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at room teinperature is too small to be of any con- 
cern. Heavy reduction did show a measurable effect 
on the temperature dependence of Young^s modulus 
at high temperature [12] but this is outside the scope 
of our present considerations. Our specimens were 
heated for about 24 hr in flowing O2 at about 800 °C 
and we consider them to be effectively stoichiometric. 

The second type of defect appears in all rutile 
single crystals examined by the authors and no 
method of removing it is known. Back-reflection 
Laue patterns were taken at random on many speci- 
mens ; perhaps one fourth of the patterns show double 
spots indicating a change of orientation of 1 or 2 deg. 
If this were a systematic, progressive change of 
orientation along a rod, one end would differ greatly 
in orientation from the other and the crystals would 
be useless for our work. Fortunately, as table 1 
shows, this is not so. Laue patterns were taken at 
the center and at both ends of each rod. Duplicate 
patterns were taken on specimens 21 and 22 to indi- 
cate the precision of the method. The resulting 
standard deviations were 0.12° for d and 0.16° for 
where 6 and <t> are the usual spherical polar angles. 
The colatitude, 6, is the angle between the rod axis 
and [001]; the angle </> is the angle between [100] and 
the projection of the rod axis on the (001) plane. 
The orientation apparently shifts back and forth in 
a random manner from one small angle boundary to 
another. The average values of 6 and </> were used 
for each rod in the calculations. It is hoped that 
better temperature, atmosphere, and powder feed 
control during crystal growth would reduce the 
number of the boundaries and this is being attempted 
at the National Bureau of Standards and elsewhere. 

Both the Linde Co.^ and the National Lead Co. 
were asked for high-purity single crystal specimens 
and both graciously gave us some of their production. 
The National Lead specimens were in the form of 
boules intended for gem stones, and rods cut from 
these proved to be too short for our use. Some of 
the Linde rods were long enough to permit cylinders 
of useful length to be ground from them. These 
rods are listed as numbers 21 through 30 in table 1 
and our first set of elastic constants [8] was calculated 
from results on these rods. 

The manufacturer's analysis for a typical Linde 
rutile single crystal indicates the following oxides 
at each stated percentage: O.OI-B2O3, ZnO, Zr02, 
SbsOs; O.OO5-AI2O3, U2O5, re203, SrO, MoO, PbO, 
Sn02; 0.003-SiO2, CaO, Cr203, C02O3, BaO; 0.002- 
MnO, NiO; 0.0005-MgO, Ag20, CuO. The lack of 
any rods with high 6 values in this set is apparent 
from table 1. The Linde Co. then made a special 
series of growth runs and produced rod number 44 
with a very useful orientation. Additional rods 
were also needed and a successful effort was made 
by W. S. Brower and S. F. HoUey of the National 
Bureau of Standards to grow rods parallel to the 
[100] axis. The purest available Ti02 powder from 
National Lead Co. was used. The manufacturer's 
analysis of the impurities in percent is: 0.02-SiO2; 
0.01-Nb; less than 0.005-W; 0.001-Fe2O3; less than 
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Table 1. Properties of "single crystal'^ rutile specimens 



Speci- 
men 
number 



Length 



End 



Center 



All angles in degrees 





g 


cm 














21 


2. 6317 


13. 236 


/ 12.5 
\ 12.3 


2.6 
2.2 


}n. 


2.1 


r 14.3 
1 14.2 


1.8 
1.5 


22 


2. 2113 


13. 178 


/ 37.9 
1 38.0 


12.5 
12.5 


} 38.9 


10.2 


/ 40.1 
I 40.2 


8.0 

8.1 


23 


2. 8205 


12. 168 


14.0 


7.3 


16.1 


7.0 


14.8 


6.3 


24 


2. 3654 


11. 714 


43.3 


34.1 


45.6 


37.9 


44.7 


33. 6 


25 


2. 3770 


11. 178 


11.2 


37.9 


11.4 


29.0 


11.9 


36.3 


26 


1. 2718 


7.534 


8.0 


26.6 


7.8 


22.8 


8.4 


26.1 


27 


1. 2769 


6.986 


14.6 


18.8 


14.2 


19.0 


15.0 


20.9 


28 


1. 2561 


6.080 


39.0 


40.0 


42.4 


40.0 


43.0 


45.0 


29 


0. 8562 


6.004 


12.9 


7.8 


14.0 


5.6 


13.7 


6.0 


30 


1. 1523 


5.816 


58.7 


24.0 


58.5 


24.2 


59.3 


23.7 


44 


15. 0251 


9.848 


88.8 


44.1 


88.4 


43.4 


90.0 


44.5 



NBS rods 



41 


8. 0767 


8.508 


12.9 


9.2 


12.9 


8.0 


13.5 


9.3 


42 


5. 8920 


6.256 


7.0 


20.0 


6.9 


21.0 


6.2 


19.0 


49 


7. 1208 


11. 128 


87.8 


0.0 


87.8 


0.0 


88.0 


1.5 


50 


21. 5907 


17. 107 


88.8 


5.4 


89.5 


5.4 


90.0 


5.0 


51 


16. 0263 


17. 634 


86.2 


0.4 


86.2 


1.0 


86.8 


0.9 



O.OOI-AI2O3, Sb203, Pb, V, Cr; less than 0.0005-Hg, 
less than 0.0003-Cu; less than 0.00005-Mn. Spec- 
troscopic analysis of single crystals after growth 
showed no significant change in the major impurities. 

3. Relation of Young's Modulus and Shear 
Modulus to Elastic Compliances and to 
Orientation 

In this section the elastic moduli measured in 
resonance experiments are expressed in terms of 
the elastic compliances and the orientation. Tlie 
method of solving these equations for the compli- 
ances is then described. The theory and experiment 
closely parallel work on corundum [11] and therefore 
only a brief description will be given. 

In relating the elastic compliances to Young's 
modulus, it is necessary to distinguish between 
the ^'free'' Young's modulus, Yf, and the ''pure" 
Young's modulus, Yp, when analyzing flexural or 
torsional tests. The free Young's modulus is the 
value obtained when the specimen is completely 
free to deform elastically under the applied tensile 
stress. The pure Young's modulus is the value 
obtained when the specimen is tested in flexure and 
is prevented from twisting. In an isotropic medium 
the free and pure moduh are identical. Calculation 
of Young's modulus from fiexural vibrations of 
slender, cylindrical rods of nonisotropic material 
corresponds to measurements of Yf. For such rods 
the modulus determined from torsional resonance is 
the pure shear modulus, Gp. The validity of these 
results is discussed by Brown [13]. 

The following treatment develops the relation be- 
tween the measured quantities, Yf and Gp, and the 
elastic compliances, Stj, expressed in the standard 
rectangular coordinate system for tetragonal crys- 
tals, 3:1X2X3, described by Nye [14]. It is convenient 
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to introduce a rotated coordinate system, a^Jxa-r.}, with 
the x-s axis alono- the rod axis. The free ehxstic niod- 

Bhited 
hnate 



uli for measurements aloni^* tlie rod axis are rel 
to the compliances expressed in tlie primed coord 
system by 



and 



2 '~Gf 



0) 



(2) 



Brown [13] and Hearmon [15] i^ive the following 
equation relating the pure shear inochdus, dp, to 
tlie free shear modulus, (i/. 



where 



/ ^34 ~r '^35 

^33 ("^44 + ^55) 



(3) 
(4) 



Equations (2), (3), and (4) can be coml)iiio<l lo give 



1_1 I {s^i+s',l)Y, 
Gr ^/ 2 



(5) 



It is now necessary to express the primed (juantities 
in eqs (1), (2), and (5) in terms of the unprimed 
compliances and the orientation. This is done for 
eqs (2) and (3) by writing out the tensor transfor- 
mations and expressing the direction cosines in terms 
of the colatitude, B, and the azinmth, (/>. In other 
words, 6 is the angle between x^ and x^^ and (/> is the 
angle between Xi and the projection of x'^ on the 
Xi:r2 plane. The tensor transfoi'mation gives 

^33=-^'ii sin^ ^+.s'33 cos^ ^+(2.s']3 + '^'44) ^hr^ 6 cos^ 6 

+\ (^'i2-.^ii+y) sin^ sin*^ 20, (()) 

and 

?^^^^=l s,, sin^ d+\ ,944(1 +cos2 e)+2{sn+s,, 

— S44— 2si3) sin2 6 cos^ ^+(^s„ — .^12— ^jsin^ B sin^ 2</). 

(7) 

The directions of the x[ and X2 axes have not been 
specified except for the requirement that x[x^x^ form 
a right handed rectangular system because eqs (1) 
and (2) are independent of their directions. This 
independence also holds for (5) because ^34+535 de- 
pends only on the (hrcction of x'^ as Hearmon [15] 
states. It is convenient to write .934 and .s^^s sep- 
arately, however, and these quantities are not spec- 
ified by the (hrection of x'^ alone. For simplicity 
X2 was chosen to lie in the XiXo ])laiie and x[ is tiien 
uniquely specified. Tlie tensor transformations give 



"^34 






«11 j s 



«ii ) siir' 6 sill 4</), 



(8) 



and 

•S'35 = 2(.S'h+'^33 — ^44-2-S^13) slll^ B COS ^+(2.9i3 + .94i 

-2,^33) sin B cos B^(^s,,-\-'^-s,^ 

Xsin^ (9 cos 6> sin2 2(/). (9) 

A set of elastic compliances can be determined 
from the preceding equations in the following man- 
ner: Assume that values of Yf and G^ have been 
determined on at least four rods. Equation (6) is 
written for each value of Yf and the resulting system 
of simultaneous equations is solved for -Sn, ,933, 
2'S'i3+<5^4, a,nd 6'i2 — -s'i 1+^44/2. These results are used 
in eqs (8), (9), and (7) to calculate values of Of 
from the measured values of G^. Equation (7) is 
written for each value of Gf and the resulting system 
of simultaneous equations is solved for .S44, 



'^'11 r'^33" 



-2,9,3, and 



'^'11 '^*12 '^'06, 



/2. In this way 



two independent values of .s'12 and ,s'i3 are obtained 
which may be compared for consistency. If more 
than four specimens are used, a good check on the 
consistency of the whole set of ,s'/; values is then 
available. 

4. Relation of Elastic Moduli to Resonance 
Frequencies 

The values of Yf and G ,, needed for the calculation 
of the elastic compliances were obtained by resonance 
frequency measurements on slender, cylindrical rods. 
Young's modulus was calculated from the longi- 
tudinal resonaTice frecjuency^ for the five longest rods 
and from the flexural resonance frequency for all 
of the 16 rods used hi this investigation. The shear 
modulus was calculated from the torsional j'csonance 
frecjucMicv of all IG rods. For both of the Young's 
modulus calculations, the ec^uation relating the 
resonance frequency to the appropriate modulus is 
approximate but the approximations are very goocL 

For longitudhial vibrations. Young's modulus 
of an isotropic nie(hum can be calculated from 
Rayleigh's equation [16] which can be written 



F=4p/ 



(2^2\ 2 



(10) 



where p is the density in g/cni^, / is the length in 
cm, / is the fundemental longitudinal resonance 
frequency in cycles per second, a is Poisson's ratio, 
r is the radius in cm and Y is Young's modulus in 
dynes/cm^. The term in parentheses is tlie Rayleigh 
correction term for the finite thickness of the rod 
and neglects higher powers of /'//. For the values 
of rjl used in tlie present work the difference be- 
tween values given by Rayleigli's equation and a 
more accurate treatment by Bancroft [17] is nuich 
less than the experimental accuracy. The correction 
term should be modified for noniso tropic material, 
but a good estimate of the correction can be obtained 
by using an average value of Poisson's ratio in the 
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present equation. Taking (7=0.30 and r//=:0.0113 
for the shortest of the five rods used in longitudinal 
resonance gives a Rayleigh correction factor of only 
1.000006. ^The values of Y calculated from (10) 
should be very accurate despite the fact that rutile 
is noniso tropic. 

For flexural vibrations of a cylindrical rod of 
isotropic material the best existing theory seems to 
be Pickett^s [18] differential equation.^ The result 
obtained can be expressed as the result for a rod 
of infinitesimal thickness multiplied by a dimen- 
sionless correction factor, T. This result is 



r=0.31547p "-^ T 



ly^ 



(11) 



where the symbols have the same meaning as in 
eq (10) except that here / is the flexural resonance 
frequency. Tefft [22] has calculated a table of T 
from Pickett's differential equation as a function of 
rjl and Poisson's ratio. Fortunately the T values 
are nearly equal to one and depend very little on 
the Poisson's ratio value for small rjL The values 
of T for (7=0.30 were used for all calculations. The 
whole subject of the determination of elastic moduli 
of isotropic materials has been summarized [23] 
recently. 

For torsional vibrations of an isotropic cylindrical 
rod, the equation 

G=^^pr~f (12) 




Figure 1. Method used to measure torsional resonance 
frequency of short rods (about 6 cm long) . 

The fine phosphor bronze springs have sufficient flexibility to permit relatively 
free specimen vibration and have sufficiently low damping not to cause excessive 
power loss. 



2 This equation should not be confused with an approximate differential equa- 
tion derived by Timoshcnko [19] and studied by Goens [20] and by Pickett [21]. 



is rigorously true. Here / is the torsional resonance 
frequency and the other symbols have their previous 
meanings. The significance of this equation for 
crystalline cylinders has been considered by Brown 
[13] and by Hearmon [15]. The details are compli- 
cated, but the result is, as previously stated, that 
for sufficiently slender rods eq (12) gives G^ and 
eqs (10) and (11) give Yf. 

The density value of 4.250 g/cm^ reported by 
Swanson and Tatge [24] was used. 

5. Method of Determining Resonance Fre- 
quencies 

The general procedure for measuring resonance 
frequencies of slender cylinders has been described 
previously [11, 23]. The specimen is suspended by 
fine threads tied in such a manner that free vibration 
is unimpeded; i.e., the threads are tied near the nodes 
of the vibrational mode being investigated. Vibra- 
tion is excited either by air drive with a loudspeaker 
or by driving one of the suspending threads with a 
transducer. The resonance is detected by a pickup 
attached to another thread. The resonance fre- 
quency is measured with a crystal controlled counter 
having an accuracy of ±0.1 c/s. 

This method was successfully used in the present 
work for all resonance frequencies except the torsional 
frequencies on rods 26 through 30. In the case of 
torsional vibration the resonance frequencies are so 
large and the amplitude of motion so small that the 
resonances could not be detected using cotton or 
silk thread. A special procedure, shown in figure 1, 
was used. Fine springs, made from phosphor bronze 
wire, were used to couple the driver and pickup to 
the specimen. These springs were flexible enough 
to permit free vibration of the specimen, but had 
sufficiently low damping to permit adequate power 
transfer from the driver to the specimen and from 
the specimen to the pickup to allow resonance to be 
excited and detected. In this way the torsional 
resonance frequencies of specimens 26 through 30 
were measured. 

The value of Young's modulus calculated from the 
longitudinal resonance frequency was considered to 
be slightly more accurate than the value calculated 
from the flexural resonance frequency and was used in 
subsequent calculations when it was available (i.e., 
for those rods for which the longitudinal resonance 
frequency could be determined). 

6. Results and Discussion 

The flexural and longitudinal resonance frequencies 
are given in table 2 together with the reciprocal 
Young's modulus values calculated from these 
frequencies. Table 3 gives the corresponding tor- 
sional resonance frequencies and reciprocal shear 
modulus values. The values of 1/F/, 6, and ^ were 
used to write eq (6) for each rod. The resulting set 
of 16 simultaneous equations in four unknowns was 
then solved for the four linear combinations of Sij 
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Table 2. 


Young's 


modulus calculations 






Resonance 
flexural 


Frequency 
longitudinal 


11 Vf 




Specimen 
number 


Flexural 


Longitudinal 


ir. 




10-13 cmVdyne 



Linde rods 





c[s 


els 








21 


1166. 5 


35544 


2.655 


2.658 


-0. 008 


22 


943.0 


31271 


3.490 


3.464 


-.014 


23 


1485. 2 


38514 


2.671 


2.678 


-.012 


24 


1393. 


37325 


3.080 


3.077 


+. 030 


25 


1703. 8 


42399 


2.614 


2.619 


-.023 


26 


3338. 2 
3980. 8 




2.614 
2.689 




-.003 


27 






+.011 


28 


5317. 5 




2.961 




+.024 


29 


4748. 3 




2.702 




+.032 


30 


5001. 8 
5693. 




3.831 
2. 7295 




-.050 


44 




2.699 


+.005 









41 


6107. 6 
11246 
1840. 6 
1111.6 

882.3 




2. 655 
2.617 
6. 824 
6. 617 
6.714 




-0. 018 


42 






-.007 


49 
50 
51 


2(>402 
17408 
16740 


6.814 
6.632 
6. 750 


+.037 
-.018 
-.007 



Tri,= (l/F/)observe.i—(l/y/)caicuiated where Yf computed from longitudinal reso- 
nance was used when available. 



Table 8. Shear modulus calcidations 



Specimen 
number 



Tor.sional 
resonance 
frequency 



1/G, 



VGf 



10-13 cmVdyne 



Linde rods 





els 








21 


20206 


8.213 


8. 236 


-0.046 


22 


19496 


8. 885 


9. 301 


+.245 


23 


21936 


8.250 


8.281 


-.059 


24 


20370 


10. 304 


10. 447 


-.117 


25 


24070 


8. 099 


8.110 


-.138 


26 


35737 


8. 115 


8.121 


-.037 


27 


37942 


8. 372 


8. 396 


+. 057 


28 


39340 


10. 282 


10. 298 


-.117 


29 


43920 


8. 459 


8.483 


+.188 


30 


42640 


9. 580 


10. 560 


-.014 


44 


20260 


14. 778 


14. 796 


+.047 



41 


31319 


8.285 


8.307 


+0. 023 


42 


43016 


8.123 


8.127 


-.004 


49 


26604 


6.712 


6.718 


+.019 


50 


17215 


6.782 


6.947 


-.011 


51 


16844 


6.668 


6.686 


-. 035 



TF,= (l/(?/)ob.e 



-(l/G/)caIculated. 



appearing as coefficients in eq (6). The process of 
solving an ovcrdeterniined set of simultaneous linear 
equations foi* the least squares best estimates of the 
coefficients and their standard deviations is a 
straight foi'ward calculation and is described, for 
example, by Scheffe [25]. This calculation was done 
on an automatic computer and the resulting values 
of the coefficients are given in the first half of table 4. 
These coefficients were then used in eq (6) to calcu- 
late a value of 1/F/for each rod. The deviations of 
the calculated from the observed values of 1/F/, 



labeled Wy, are listed in table 2 and provide an 
indication of how well eq (6) fits the data. These 
r(^sults were then used in eqs (8), (9), and (5) to 
o})tain IjGf from IjG-p for each rod. The values of 
IjGf were tluni used in eq (7) and a least squares 
solution for tlie coefficients was carried out as was 
(lone for eq (6). TJie resulting values of the coeffi- 
cients are given in the second \\M of table 4. Devia- 
tions, Wg, between calculated and observed values 
of IjOf are listed in table 3. 

Table 4. Least sauares best estimates of parameters in elastic 
moduli equations 

From equation for 1/F/ 



sn 


6.788±0.015X10-i^ cm^dync 


S33 


2.592± .011 


2513+S4-1 


6.466± .062 


sm 


8.200± .059 


Si\ — sn — - 





From equation for 1/(7/ 



*'ll+«33-'Sl4— 2«13 
Sll— «12— W2 



8.072it0.048 
5.302± .140 
2.880± .157 
8.077± .132 



The deviations are standard deviations of the coofficionts obtained from least 
squares calculations. 

The results presented in table 4 give Su, %3, '*<44, and 
Sqc directly. A valuc^ of .Si2 and -s'l^ can })e calculated 
from tlieir appearance in tlie 1/)^ coefficients. A 
second valu<' of eacJi can be calculated using the 
IjGf coeflicicMit in wJiich thev appear. For Sn the 
two values are — 4.()63±0.088 and -:].940±0.1i:^; 
for Su tlie result is -0.803 ± 0.0:^9 and - ().78() ± 0.069 
The results are thus self consistent. These values 
were weighted according to the reciprocals of their 
variances and tlie values of ^12=— 4. 017 ±0.069 and 
.s'i3= — 0.799 ±0.034 were obtained. 

The Cij values can be calculatcMl by inversion of the 
matrix of s^ values. The equations for passing from 
dj to Sij are given by Nye [1/f] and it is easy to show 
that the equations for obtaining Cij from -s*,-^ in this 
case are obtained by interchanging ^'s and 6*'s. The 
result is 

1 

^11 = 



1 r^3 , L_1 



^33 = 



C44- 



^66 = 



s 

6*44 
1 



s^ 



^^ 2Ls 6'ii— 512J 



Cn=- 



^13 
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where 



^ — "^S.sC'^ll I '^12) ^^13- 



The calculations were made b}" substituting for the 
Sij the directly computed quantities of table 4. For 
each of the constants Cn, C33, C12, and C13 the two 
values were obtained and a weighted mean computed 
as in the case of ^12 and ^is. 

The Sij and Cfj values are listed in table 5 for com- 
parison with the work of other investigators. The 
present Ctj show no statistically significant differences 
from the values calculated by Birch. Using the 
simple test of twice the standard deviation we find, 
however, that the ^n and ^12 values of Birch are 
significantly different from our results. This is 
probably caused by the fact that the equations for 
^11 and S12 (which have the same form as the above 
eqs for Cu and c^) both involve the difference Cn — c^. 
There is thus some reason to prefer the present ^n 
and ^12 values as they are more directly determined; 
the same argument favors the Cu and C12 values of 
Verma-Birch. Our linear and volume compressibil- 
ity values are also given in table 5 and are consistent 
with Bridginan and Verma-Birch. 

It is interesting to compare the elastic properties of 
rutile single crystals with those of polycrystalline 
rutile. Huntington [26] discusses the problem of 
calculating elastic moduli of a pore-free polycrystal- 
line solid from the single crystal elastic constants. 
Results calculated from theories of Reuss (which 
gives a lower bound) and of Voigt (which gives an 
upper bound) are given in table 6. The bulk 
modulus is determined within narrow limits but the 
Young's modulus and shear modulus values are not 
so well determined. This is a reflection of the fact 
that, as shown by figures 2 and 3, these two moduli 
depend strongly on orientation. 

Table 5. Elastic parameters of rutile 



Cii in 1012 dyne/cm2 


Sij in 10-13 cm2/dyne 




Bridg- 
man » 


Dayal 

and 
Appala- 
narasim- 

ham b 


Vick and 
Hollander " 


Verma 
recalcu- 
lated by 

Birch d 


Present 
work e 


Cll 
C33 
C44 

C66 

Cn 
Cn 
sn 

533 
544 
S66 
512 

Sl3 


"" i.'89" 
1.04 

4.82 


3.005 

1.9 

1.324 

1.761 

1.76 

1.36 

5.8 

8.9 

7.6 

5.7 
-2.2 
-2.5 


2. 48±. 08 
4. 52±. 08 
1. 20±. 03 
1.6 ±.1 
2.0 ±.1 
1.4 ±.1 
11.8 
2.7 
8.3 
6.2 

-9.0 

-0.86 


2.73 
4.84 
1.25 
1.94 
1.76 
1.49 
6.55 
2.59 
8.00 
5.16 
-3.76 
-0.86 
1.93 
0.87 
4.73 


2.660±.066 

4. 699±. 081 
1.239±.007 
1. 886±. 050 
1.733±.071 
1.362±.081 
6.788±.015 
2.592±.011 
8. 072±. 048 

5. 302±. 140 
-4.017±.069 
-0. 799±. 034 

1 965_j_ 069 


S33+2si3 






0. 994±. 067 


2S] i+5m+2si 2-l-4i!''i ■? 






4. 911± 166 











a Measured statically. Present adiabatic values calculated from Bridgmen's 
isothermal values. 

b Cii computed from spectroscopic observations, sa computed by matrix 
inversion. 

c Diagonal Ci,- computed from pulse velocity measurements. C12 and cn com- 
puted using Bridgman's linear compressibilities, sa by matrix inversion. 

d Cji computed from pulse velocity measurements, sa by matrix inversion. 

e Sij computed from resonance frequency measurements, ca by matrix invercion 
as explained in the text. The values for the linear and volume compressibilities 
given in the last three rows are not exactly the same as would be obtained by direct 
computation from the Sij above, because the weighting is different for a combina- 
tion of constants than for a single constant. The deviations shown for the present 
values are standard deviations for the compliances and constants obtained from 
least squares calculations. 



Table 6. Elastic moduli for polycrystalline rutile computed 
from Reuss and Voight theories * 



Young's modulus 
Shear modulus 
Bulk modulus 



Polycrystalline 



Reuss 
theory 



Voight 
theory 



Single crystal 



All units 1012 dyne/cm2 



2.555 
0.990 
2.025 



3.116 
1.244 
2.103 



Orientation dependent see figure 2. 
Orientation dependent see figure 3. 
2. 025. 



a Using values of the Sij and cij given in the last column of table 5. 

The writers have previously used the resonance 
method to determine the six elastic compliances, 
Sij, of corundum [11] with standard deviations of 
about 0.1 percent for the diagonal compliances and 
about 1.0 percent for the off -diagonal compliances. 
The present uncertainties in the diagonal constants 
range from 0.2 percent for Sn to 2.6 percent in 6^66- 
For the off-diagonal constants the values are 1.7 
percent for ^12 and 4.3 percent for 813. These results 
are somewhat disappointing when compared with 
the precision achieved with corundum. Rutile 
difl^ers from corundum in being more anisotropic 
but this should not cause such an increase in standard 
deviation if the orientation were constant through- 
out the specimen. The principal difficulty probably 




Figure 2. Young^s modulus, Yf, as a function of orientation, 
calculated from eqs (1) and (6), using the least squares best 
estimates of the coefficients given in table 4- 
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Figure 3. Shear modulus, Gf, as a function of orientation, 
calculated from eqs {2) and (7), using the leak squares best 
estimates of the coefficients given in table 4- 

is assoeicated with the small angle l)()uii(laries. If 
speeinuMis free of these boundaries arc* ever produced 
it would he interesting to repeat both the pulse 
velocity and resonance measurements and attempt 
to obtain agreement to 0.1 percent in the diagonal 
constants. It is probably not worthwhile to attempt 
to improve the agreement using existing specimens. 
The present results are accurate enough for most 
purposes and the general agreement of different 
observers leaves no question of serious error. The 
value of C33 calculated by Dayal and Appalanara- 
simham from spectroscopic data is in serious error, 
as previously noted [9, 10], and an extension of their 
calculation to correct this value might be worthwhile. 

Note. — -The writers were recently informed by Oilman 
[27] of unpublished work by himself and B. Chick on the 
elastic constants of rutile. They obtained C33=4.75X lO^'^ 
dynes/cm^ and C44 = 1.23X 10^^ dynes/cm^. Both values are 
in good agreement with the results of the present work. 
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